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Abstract. The long-term variability of the multiwavelength blazar emission can be interpreted in terms of orien- 
tation variations of a helical, inhomogeneous, non-thermally emitting jet, possibly caused by the orbital motion 
of the parent black hole in a binary system (Villata & Raiteri 1999,). The helical-jet model is here applied to 
explain the quasi-periodic radio-optical light curves and the broad-band spectral energy distributions (SEDs) of 
the BL Lac object AO 0235-1-16. Through a suitable choice of the model parameters, the helix rotation can well 
account for the periodicity of the main radio and optical outbursts and for the corresponding SED variability, 
while the interspersed minor radio events could be interpreted as due either to some local distortions of the helical 
structure or to other phenomena contributing to the source emission. In particular, the probable existence of flow 
instabilities provides a viable interpretation for the non-periodic features. 

Key words, galaxies: active - BL Lacertae objects: general - BL Lacertae objects: individual: AO 0235-1-16 - 
galaxies: jets - galaxies: nuclei - quasars: general 



1. Introduction 

Blazars, namely BL Lacertae objects and flat-spectrum 
k> , radio quasars, belong to the class of active galactic nuclei 
^ ■ (AGNs). According to the unified model of AGNs (see 
5^ ! G-g- Urry & Padovani 11995(1 . their emission is dominated 
by radiation produced in relativistic plasma jets, which 
are oriented at small angles with respect to the line of 
sight; the radiation is therefore strongly beamed towards 
the observer. 

Details of the jet origin are still unknown, but their 
formation is thought to be triggered by the presence of a 
supermassive black hole surrounded by an accretion disc, 
possibly belonging to a binary black hole system (BBHS) 
hidden in the centre of the AGN. On one hand, this sce- 
nario is suggested by the elliptical morphology of BL Lac 
host galaxies, which are believed to originate from merg- 
ing phenomena between spirals, naturally leading to the 
formation of massive binary black holes (Begelman et al. 
[TUm Wilson & Colbert [T^ Huges & Blandford l^nUSjl . 
On the other hand, several observational evidences such 
as bending, misalignment, wiggling and precession of jets 
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(Begelman et al. 11980? Camenzind & Krockenberger[ 
Kaastra & Roos il992. Conway & Wrobel . 19951 Villata et 
al. 119981 Villata & Raiteri 119991 and references therein; 
Abraham '2000': Britzen et al. l2001|l . often associated with 
knots superluminally moving along different-scale curved 
trajectories, as well as the periodicity discovered in the 
multiwavelength light curves of some of these sources, have 
in some cases been interpreted in terms of helical struc- 
tures tightly related to BBHSs. 

The helical-jet model proposed by Villata & Raiteri 
{1999| describes how orientation variations of the 
different-frequency emitting parts of a curved jet with re- 
spect to the line of sight, possibly caused by the orbital 
motion of the parent black hole in a BBHS, can be the 
cause of the observed changes in the spectral energy dis- 
tribution (SED) of a blazar. The emission from the jet 
is non-thermal: the relativistic electron population is re- 
sponsible for producing both the low-energy (from radio 
to UV-X-rays) synchrotron radiation and the high-energy 
(up to 7-rays) one through inverse-Compton (IC) scat- 
tering of the same synchrotron photons (synchrotron-self- 
Compton process, hereafter SSC). 
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This model has already provided an interpretation for 
the huge X-ray spectral brightening of Mkn 501 (Villata 
& Raiteri J199S|I, the low-energy SED variations of S4 
0954-f 65 (Raiteri et al. ll999|l . and the changes in the over- 
all SEDs of S5 0716-^71 (Ostorero et al. lMTTll and ON 231 
(Sobrito et al. l2(J(Jl|l . The modelling of a constant helix ro- 
tation allows us now to simulate the long-term behaviour 
of the multifrequency light curves together with the SED 
time evolution. If the data sampling is good enough and 
the source emission exhibits a periodic (or quasi-periodic) 
behaviour, one can apply the model and find the set of 
parameters which best reproduces the observed multifre- 
quency data, thus providing both a test for the model and 
a description of the corresponding jet features. 

In this paper the model, fully described in Sect. 2, is 
apphed to the case of the BL Lac object AO 0235-1-16, 
whose radio (and optical) light curves have recently re- 
vealed a ^ 5.7 year quasi-periodicity (Raiteri et al. '2001). 
Model light curves and SEDs are compared with observa- 
tional data in Sect. 3. Possible interpretations of minor, 
non-periodic events are presented in Sect. 4. A final dis- 
cussion is performed in Sect. 5. 



The jet is assumed to be inhomogeneous, similarly to 
that of the model of Ghisellini & Maraschi 1)19891 see also 
Ghisellini et al. HMSl Maraschi et al. HW^- each slice 
of it emits synchrotron radiation in a range of frequen- 
cies between u'^ and v!. ^^^^ (primed quantities refer to 
the plasma rest reference frame), which are supposed to 

decrease as power laws as the distance I from the apex 
1. 



mcreases 



c(0 



.i„(0) ( 1 + ^ 



c(0) 1 



ci > 0, 



C2 > 0, 



(4) 



(5) 



where v'^ niin(O) ^^'^ ^'s max(O) ^'^^ the values at ^ = z = 0, 
and I2 are length scales. In the following we assume that 
^'s,min(0) = ^'s,max(0) = Kfi- The jet cmits also higher- 
energy radiation, via IC scattering of soft synchrotron 
photons off high-energy electrons (SSC process), in a range 
of frequencies 



^c,min(0 ^Tr^iTi^* 



mm s, mill 



(0 



(6) 



2. The helical-jet model 

The helical-jet model by Villata & Raiteri H1999fl foresees 
that, in a BBHS scenario, the emitting jet is bent because 
of the orbital motion and by the interaction of its magne- 
tized plasma with the surrounding medium, twisting in a 
rotating helical-shaped structure. 

Since in Villata & Raiteri H1999|l only the synchrotron 
emission was taken into account, we recall here the main 
formulae, adding those expressions relevant to the high- 
energy SSC emission (see also Raiteri et al. I2004|l . 

Let us consider a helical-shaped jet with the helix axis 
taken as z-axis of a 3-D reference frame. The helix has 
a pitch angle C and the line of sight forms an angle ij) 
with the z-axis. The length of the helical jet path can be 
expressed as 



l{z 



< z < 1 , 



cosC ' 

and covers an azimuthal angle 
ip{z) = az, 



(2) 



where a is the total angle covered by the helix. The jet 
viewing angle at z is 



cos 6{z) — cos -0 cos ( + sin ip sin ( cos( 



(3) 



where (p is the "rotation angle" , i.e. the azimuthal differ- 
ence between the line of sight and the initial direction of 
the helical path. 

The jet geometry described above is the apparent 
shape of the plasma flow; the different light travel times of 
each jet region and the plasma velocity components due 
to the rotation of the helical structure have been taken 
into account. 



.(0 



c 



(7) 



where the electron Lorentz factor 7 is supposed to vary 
between 7inin = 1 and 



7max(0 = 70 1 + 



, > , 



(8) 



being 70 = 7max(0). 

Actually, in order to take the Klein-Nishina effect into 
account, an upper limit is set to v'^ max(Oi the form 

f 4 m c'^ 1 

'^c,max(0 = 7max(0 | 3 7max (O^^lmaxCO : —J^ > ■ (9) 

The observed flux density at frequency v is assumed 
to be 



(1) F^{iy) oc 5^1^-°"' , 



(10) 



where ao is the power-law index of the local synchrotron 
spectrum (fixed to be ag — 0.5) and 



[r(l -/3cos6i)]^i 



(11) 



is the beaming or Doppler factor, /3 being the bulk velocity 
of the emitting plasma in units of the speed of light, T — 
(1 — /3^)~^/^ the corresponding bulk Lorentz factor, and 9 
the angle between the velocity vector and the line of sight, 

^ In the Ghisellini & Maraschi (119891 model, what is 
parametrized and described with power laws are the physi- 
cal quantities of the jet, on which the emitted frequencies de- 
pend. Since we are interested in Doppler factor variation ef- 
fects (whatever the physical details can be), we prefer a more 
direct, phenomenological approach, in order to lower the pa- 
rameter number. If needed, one can go back to the physical 
characteristics. 
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which varies along the hehcal path as shown in Eq. 
hence causing a dependence of S on the z coordinate. 

Notice that the local viewing angle 9 of the jet reaches 
a minimum when (p — az is zero or a multiple of 360°, 
implying the maximization of the beaming effect for fre- 
quencies emitted from the corresponding jet segments. 

By introducing an intrinsic flux-density dependence on 
I, the observed synchrotron and IC flux densities at fre- 
quency coming from a jet slice of thickness d^ can be 
expressed as 



dF,,e(i/) = K, 



I 



d^{l)iy-°'°dl, Cs>0, (12) 
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where and Ic are again length scales, and the constants 
Kg and Kc are considered independent of time, meaning 
that intrinsic variations of the flux are not allowed. The 
respective total flux densities at frequency v coming from 
the whole jet are obtained by integrating along all the 
portions contributing to that observed frequency: 



AzB.i(i^) 
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where Al/Az = 1/ cosC, and Azs,c:i('^) is the z interval 
corresponding to the z-th jet segment emitting radiation 
at frequencies observed as v along the line of sight, i.e. 
where the condition 



(16) 



is satisfied. An explanation of the logarithmic term in Eqs. 
((T^ and 1^3) can be found in GhiscUini & Maraschi 1989 

The total flux density at frequency v is finally obtained 
by summing the synchrotron and IC contributions. 

If the helix rotates, e.g. because of the orbital motion 
of the parent black hole in a BBHS, the ip angle, and hence 
9 [Eq. ||2Jl], vary with time; this consequently makes the 
Doppler factor 6 [Eq. (|ll|l ] of each emitting portion Al of 
the jet time-dependent, thus causing a temporal evolution 
of the observed flux density F^{v), even if the jet emission 
has not intrinsically changed at all. 



3. Modelling the periodic long-term behaviour of 
AO 0235+16 

The recently discovered periodic behaviour of radio (and 
optical) fight curves of AO 0235-1-16 (Raiteri et al. lW)T)l 
and the wealth of historical spectral coverage of this source 
suggested us the advisability of applying the rotating 
helical-jet model, intrinsically periodic, to model the long- 
term trend of its multifrequency light curves and to de- 
scribe the corresponding broad-band SED time evolution 
(preliminary results of this work have already been pub- 
lished in Ostorero et al. l2003c3 and Ostorero et al. l2003bjl . 

With this purpose, we assembled the historical light 
curves and composed the multi-epoch SED of AO 
0235-1-16, by collecting data at all available frequencies 
from the literature and homogenizing them through suit- 
able corrections for extinction effects. In particular, opti- 
cal and infrared magnitudes were converted into fluxes by 
using the absolute calibrations of Bessel H1979|l for stan- 
dard UBVKF and Cousins' RI filters , Wam steker itTMT^ 
for standard JHL bands, and Allen (jl973|l for standard 
RIN filters; infrared-to- ultraviolet dereddening was com- 
puted by adopting the extinction laws of Rieke & Lebofsky 
fl985') and Cardelli et al. fl9891, assuming a Galactic ex- 
tinction Ab — 0.341 mag (from NED). In the X-ray band, 
absorbed fluxes were corrected for total absorption fading 
by using the XSPEC procedure. 

3.1. The radio-optical light curves 

The most time-extended and best-sampled light curves of 
the source are the optical ones (see Fig. la), and the ra- 
dio ones at 14.5 GHz (Fig. 16), 8.0 GHz (Fig. Ic) and 
4.8 GHz (Fig. Id) from the University of Michigan Radio 
Astronomy Observatory (UMRAO). All these curves ex- 
hibit strong variability: many large-amplitude outbursts 
lasting from several months to a few years witness an in- 
tense source activity. From the long-term time analysis 
of the best-sampled light curve at 8.0 GHz, Raiteri et al. 
I)2001|) found the signature of a quasi-periodic behaviour, 
also recognizable in the other radio bands: the major out- 
bursts repeat quasi-regularly every ~ 5.7 years, many of 
them accompanied by simultaneous or slightly preceding 
optical ones. A visual inspection of the curves in the left 
panels of Fig. 1 reveals that the optical events are how- 
ever sharper, more pronounced, and sometimes more fre- 
quent than the radio ones. Indeed, Takalo et al. ( 1998) in- 
voked the existence of two different optical emission mech- 
anisms: one also responsible for the radio events (likely 
synchrotron processes in the relativistic jet), and the other 
not (probable microlensing effect of a foreground galaxy 
on radiation coming from either the accretion disc or other 
jet regions). Moreover, the radio outbursts have decreasing 
amplitude with increasing wavelength, and flux variations 
at the lower frequencies lag those at the higher ones (see 
the cross-correlation analysis performed by Raiteri et al. 
I2()()l|l : other minor outbursts are present in the radio light 
curves beside the periodic events. Finally, all outbursts 
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Fig. 1. Left panels: Observed optical (mJy) and radio (Jy) light curves of AO 0235+16. a: R optical band {R band data, 
and B band data before JD = 2449000 converted into R ones by adopting a mean colour index < B — R >= 1.65 + 0.16; 
see Raiteri et al. IMiTl for details); b,c,d: 14.5, 8.0, 4.8 GHz radio bands (UMRAO); e: 1.4 GHz (dots; Green Bank 
telescope) and 1.5 GHz (crosses; NRAO) radio bands; data have been drawn together to better define the light curve 
features; /: 880 MHz radio band (Green Bank telescope); g: 430 MHz radio band (Arecibo telescope). These light 
curves are based on data from: Owen et al. H1978|l . Owen et al. H1980|l . Ulvestad et al. 1)1981(1 . Perley H1982|l . Briggs 
ifTW). Landau et al. ^E^, Ulvestad et al. ^E^, Altschuler et al. ifBIMIl . Ulvestad & Johnston ifTMll . AUer et 
al. U2S5), Rudnick et al. ifTHH^ . Brown et al. lfT??k??| . Salgado et al. ^M^, Raiteri et al. CTIT)! and references 
therein. Right panels: Helical-jet modelling of the typical multifrequency outbursts shown in the left panels (see the 
text for a detailed description of the model parameters relevant to curves a'-g'): as the wavelength increases, the 
outbursts appear to be more and more delayed and to have smaller amplitude and longer duration, in agreement with 
observations. Grey (yellow in the electronic version) strips highlight here the time interval between the peak of the 
modelled optical outburst and the maximum of the associated lowest-energy radio event. The same strips have been 
drawn on the observed light curves of the left panels, in correspondence to the periodic outbursts, even if in more than 
one case some counterparts are poorly sampled or missing at all. 



and flaring events seem to be superimposed to a slowly 
decreasing "base level" flux (see Sect. 3.2 for details). 

Confirmations of these behaviours come from the 
lower-energy radio light curves, notwithstanding their 
poor sampling. The delayed counterparts of both the 
1982 periodic radio outburst and the 1980-81 non-periodic 
event are likely represented by the couples of peaking 
structures sequentially appearing in the curves at 1.4 GHz 
(Fig. le; dots) and 880 MHz (Fig. 1/) from the Green 
Bank telescope, at 1.5 GHz (Fig. le; crosses) from the 
NRAO, and at 430 MHz (Fig. Ig) from the Arecibo tele- 
scope. Furthermore, the outburst peak flux seems to de- 
crease with increasing wavelength more rapidly in the case 
of the periodic events than in the case of the non-periodic 
ones. Thus, while at high radio frequencies the periodic 
outbursts dominate the non-periodic ones, at low energies 
the vice versa seems to occur; this behaviour implies dif- 
ferent radio broad-band spectral indices for the periodic 
and non-periodic components. Moreover, the light curve at 
430 MHz (Fig. Ig) also exhibits the presence of a decreas- 
ing "base level" flux, on which the peaks corresponding 
to the periodic trend fade almost completely away. Notice 
also that the "base level" flux seems not to decrease with 
decreasing radio frequency (a similar behaviour was rec- 
ognized in the BL Lac object S5 0716+71; see Raiteri et 
aL lMll . 

As anticipated in the previous section, the rotation 
of an inhomogeneous helical jet may account for the pe- 
riodic behaviour of multiwavelength light curves. As the 
helix rotates, different-frequency emitting portions of the 
jet approach the line of sight, giving rise to outbursts ob- 
served at lower and lower energies as long as the rotation 
angle ((/>) increases: this phenomenon comes from both the 
jet inhomogeneity and the beaming effect, as it will be dis- 
cussed in the following. 

The radio-optical hght curve modelling was performed 
by assuming a period of 2069 days (see Raiteri et al. 2001j 
for the constant helix rotation, and by fixing (j) — at 
JD = 2451032.6341, corresponding to the 1998 outburst 
observed peak {R band). 



The strongest constraints to be taken into account in 
setting the model parameters are the observed outburst 
amplitudes, their durations, and the shortness of the delay 
of the radio events with respect to the optical ones, which 
strongly depend on the inhomogeneity of the first part of 
the jet. 

This inhomogeneity was modelled to meet the observa- 
tional constraints by replacing the general single power- 
law trends for the emitted frequencies of the helical-jet 
model [Eqs. Q and lO] with more complex laws. 

In particular, the power law describing the mini- 
mum synchrotron frequency emitted along the jet path, 
inin(0; Substituted by a combination of power laws 

for I < Zbrcak: 

<oji ( 1 + ) ' cij,ciji >0, (17) 

and with the tangent to logi/'^ ^^^{l) in / = ^brcak, for 
I > 'break J while the I's^maxlO relation is simply a lin- 
ear combination of power laws like that of Eq. ^: 

/ I \^''^-" 

'^s.OJI ( 1 + 1^ j , C2,I,C2,II > 0. (18) 

Both relations are represented by the continuous (blue) 
curves of Fig. 2 (the parameter values are given in Table 
1, Col. A). 

The minimum emitted frequency law is very steep at 
first, implying the emission of frequencies down to radio 
bands from the immediate vicinity of the apex of the emit- 
ting jet. This allows the observed optical and radio out- 
bursts to follow one another with short time delays, as the 
contiguous jet slices mainly producing them become more 
aligned with the line of sight during the helix rotation. 
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Fig. 2. Solid (blue) curves represent the minimum and 
maximum synchrotron frequencies emitted by the mod- 
elled jet of AO 0235+16, as a function of the helical path 
coordinate I. Note the steepness of the v'^ min(0 ^^'^ the 
first jet slice and its slowlier decreasing trend right after 
the apex region. The vertical, dotted (purple) line rep- 
resents the point /break = 0.017 where the law has been 
replaced by a linear decreasing trend. The law parameters 
are visible in Col. A of Table 1. The dashed (red) curve 
shows the decrease, along the jet, of the intrinsic fading 
factor: it reduces dramatically the emission well before the 
first tenth of the path, making the jet region far from the 
apex marginally contributing to the whole emission. Only 
the portion of the path between I ~ and / = 0.1 has been 
indeed represented here. For graphic convenience, the fad- 
ing factor has been normalized to 15. 



Then, the linear decrease of log u'^ ^^^^ makes the emis- 
sion of the observed very soft radio photons possible far 
from the jet apex, before the intrinsic flux-density depen- 
dence on I, represented by the "fading factor" (l + ///s)~'^' 
[see Eq. (11)], strongly reduces the jet emission (see the 
dashed, red curve of Fig. 2). 

The results of the helical-jet modelling of the light 
curves of Fig. 1 (a-g) are displayed in the right panels 
of Fig. 1 {a'-cf). Since the model light curves are exactly 
periodic, only snapshots including outbursts are shown for 
a visual comparison between the modelled and observed 
features. Grey (yellow in the electronic version) stripes, 
defined as the time interval (~ 242 days) between the op- 
tical theoretical peak and the lowest-frequency radio one, 
guide the eye through the periodic outbursts in different 
energy bands. Besides the periodic recurrence, also the 
typical shape of the main observed optical and radio out- 
bursts, the decrease of the peak flux with decreasing ra- 
dio frequency, and the optical-to-radio flux ratio are well 
described by the model. The quasi-simultaneity of opti- 
cal and radio outbursts also is well reproduced: radio de- 
lays start from about a couple of months in the case of 
the shortest wavelengths, consistently with the results of 
the cross-correlation analysis performed by Raiteri et al. 



(|2001|l . and gradually increase up to ^ 8 months, as one 
moves towards the lowest frequencies. 

A complete summary of the parameters describing 
the geometrical and physical properties of the modelled 
AO 0235+16 jet is given in Table 1: both sets A and B 
(see Sect. 3.2 for details) give rise to the same radio-to- 
optical light curves. The emerging picture for the jet of 
AO 0235+16 is a rotating helical-shaped structure with 
a pitch angle ( — 20° and the axis oriented at an angle 
ip — 15° with respect to the line of sight, implying for each 
jet slice a maximum alignment of 5°. The plasma flows 
along the helix with a bulk Lorentz factor F = 10, and 
emits synchrotron and IC radiation, all photons mainly 
coming from a jet region close to the apex. 

In general, the parameter values reported in Table 1 
represent almost the unique set of parameters which can 
strictly account for the multifrequency light-curve (and 
SED, see below) behaviour, due to the strong constraints 
imposed by the exceptional dataset. However, some pa- 
rameters, namely the jet geometry, F, and the IC parame- 
ters, are less constrained due to two independent reasons. 
For the former ones, what is strictly constrained by the 
data are the Doppler factor variations, so that other com- 
binations of F and angles can work as well. On the other 
hand, the high-energy part of the SED is poorly sampled 
(see next subsection) and does not allow one to identify 
a unique choice of parameters. In conclusion, the phe- 
nomcnological parameters (and the corresponding phys- 
ical quantities, see footnote 1) we have here derived to 
match the data can represent the "most reasonable" jet 
features, but not the sole solution. 

3.2. The SED time evolution 

The BL Lac object AO 0235+16 has been monitored in a 
wide range of frequencies since 1970's. Radio observations 
were performed by many antennas, starting from about 
100 MHz up to 300 GHz; infrared space and ground data 
are available from submillimetric (far-IR) down to micron 
wavelengths (near-IR); optical UBVRI bands have been 
intensively monitored by many telescopes all around the 
world, while only one ultraviolet upper limit by lUE is 
available. AO 0235+16 was also detected in the X-ray 
band by several space observatories operating within the 
range 0.5-60 keV (Einstein, EXOSAT, ROSAT, ASCA, 
RXTE) starting from 1978. The 7-ray emission was inves- 
tigated through the instruments onboard the Compton 
Gamma-Ray Observatory (CGRO): COMPTEL probed 
the low-medium energy 7-ray range of the source spec- 
trum during CGRO Cycle 4 (1994-1995), and only set up- 
per hmits to the flux in the interval 0.75-30 MeV; EGRET 
explored the higher-energy 7-ray radiation (30 MeV-20 
GeV) with six pointings between 1992 and 1997, and high- 
confldentially identified AO 0235+16 as a powerful 7-ray 
emitter, providing two detections and four upper limits. 

Figs. 3 and 4 show the source multi-epoch broad-band 
SEDs, composed by collecting literature data at all avail- 
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log V [Hz] 

Fig. 3. Multi-epoch SED of AO 0235+16. At radio, infrared, optical and ultraviolet frequencies, squares represent 
historical variation ranges when joined by vertical bars and single sporadic data otherwise; data are from: Condon & 
Jauncey ifTUT^ . Ryle et al. l|TW5|l . Altschul er fc W ardle ifWH)! . Rieke et al. (fW^ . O'Dell et al. ifWTIl . Rieke et al. 
ifWTji . O'Dell et al . (jMHill, O'Dell et al. I|1978b|l . Owen et al. fMHI), Condo n et al . ifTnTHjl . Epstein et al. lfTM7| . 
Landau et al. ifTMIjl . Owen et al. {TgSO), Weiler & Johnston (1980), Baath et al. lfTMT|l . Jones et al. (p'OSl), Altschuler 

Ennis et al. ijl^g^ . Impey et al. (jIlHS), Altschuler (,1983|,) . Briggs iP^S^ . Landau et al. 1(1^83^1 . Seielstad et 
al. ifTUSHjl . Sitko et al. l(TU5^ . Ulvestad et al. ifTUS^ . Altschuler et al. CSHl, Cotton et al. (HHHII), Cruz-Gonzales 
& Hucra ifDMIl . Gear et al. p984). Holmes et al. ri984), Impey et al. (1984), Ulvestad & Johnston l|TMHl . Gear et 
al. (11933, Moles et al. (fTU5^ . Rudnick et al. (P785() . Sitko et al. ^985,) . Brindle et al. (HHHOI), Gear et al. l(TU5^ . 
Jones et al. 1(1^5^ . Edelson l(TWjl . Haddock & Aller ((TW|l . Salonen et al. (ITW|) . Smith et al. ijTWji . Fugmann 
& Meisenheimer (ITO5F). Impey & Neugebauer Steppe et al. (fEHSI, Wardle & Roberts ((TU55|l . Webb et al. 

(UnHHIl, Brown et al. fl989;) . Vetukhnovskaya (fTM^ . Impey & Tapia (jWO), Mead et al. (fTM7|l . Simonetti & Cordes 

Gregory & Condon lfTMT|l . Sitko & Sitko (f!MT)l . Edelson et al. (iJ_992), Quirrenbach et al. |1992), Takalo et 
al. Terasranta et al. ifTM^ . Gear (1993), Krichbaum et al. l(T!M|l . Steppe et al. ifBMjl . Bloom et al. lfT?M| . 

Gear et al. (P1M| . Ghosh et al. lfT!M| . Lichtfield et al. lfT!M)l . Tornikoski et al. lfT!M| . Cotton et al. ifTMTji . Renter 
et al. ifTMTji . Shen et al. (f1W|l . Nartallo et al. (1998), Takalo et al. lfT?M|l . Salgado et al. (fTMn|l . Raiteri et al. ipHTll 
and references therein. At high-energies, all data have been drawn: X-ray data (Einstein, EXOSAT, ROSAT, ASCA, 
RXTE) are from WorraU & Wilkes l|TMn|) . Ghosh & Soundararajaperumal l(TM5j) . Madejski et al. l|TM5|) . Comastri 
et al. ifTMTj) . Webb et al. iPTOjl . and the TARTARUS database; 7-ray data (COMPTEL and EGRET) are from 
Schonfelder et al. H2000|) . Hartman et al. (1999), and Raiteri et al. (|2yUl'). The shaded (yellow) area shows the "base 
level" flux (synchrotron component: dashed line; IC component: dashed-dottcd line). Different-style (and colour) lines 
represent the SED evolution modelled with the helix rotation, i.e. with different values of the rotation angle (f) (see 
the legend in the figure) , under the hypothesis that the whole X-ray emission comes from IC scattering of synchrotron 
photons off jet electrons. Model parameters are shown in Col. A of Table 1. 



able wavelengths: the historical variation ranges of the 
emission from radio to optical bands are represented to- 
gether with all the X-ray and 7-ray data. The low-energy 
part of the SED exhibits a variability which becomes less 
pronounced moving from optical bands towards soft ra- 
dio frequencies, as it is also evident from the previously 
discussed light curves of Fig. 1. A variability range and a 



clear broad-band behaviour cannot instead be easily iden- 
tified for the high-energy spectral component, because of 
the poor sampling of the X-ray and 7-ray light curves, 
and the uncertainties in data analysis, as detailed in the 
following. 

In the X-ray energy domain, the source reveals a re- 
markable variability: the 1 keV flux spans less than one or- 
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Fig. 4. Multi-epoch SED of AO 0235+16. Historical data and "base level" flux are the same of the SED in Fig. 3. 
Helical-jet model fits have been here performed by supposing the existence of a synchrotron contribution to the soft 
X-ray emission, while the hard X-rays have still an IC origin. Model parameters are shown in Col. B of Table 1. 



der of magnitude, but the spectral slope shows an intrigu- 
ing alternance of steepening and hardening in the energy 
range 0.5-10 keV, the spectral index a [F^ oc v~°') rang- 
ing from 0.41 to 2.25. In particular, the spectrum seems 
to be steeper when the source is in a brighter state, as 
already noticed by Madejski et al. H1995fl . who analysed 
the spectral changes of the source from the 1993 ROSAT 
detections to the 1994 ASCA ones: such a behaviour, in 
disagreement with a simple synchrotron cooling model for 
the X-ray emission, could reveal that in AO 0235+16 two 
different spectral components overlap in the X-ray band, 
the drop in the intensity of the softer one uncovering the 
harder one. However, the softer component, which could 
be the tail of the synchrotron spectrum, is not so clearly 
identified, because of the large uncertainties in the X-ray 
spectral indices (see the X-ray butterflies in Figs. 3 and 
4). The main problem in the data analysis concerns in- 
deed the estimate of the soft X-ray absorption. A model 
with a power-law continuum and absorption usually fits 
well the spectrum of AO 0235+16, and the best fit is ob- 
tained by letting the hydrogen column density parame- 
ter A^H free; in this case, however, A^h is larger than the 
standard value. Alternatively, one can keep A^h fixed, ob- 
taining smaller uncertainties on the spectral indices but 
a worse fit. From spectral studies, Madejski et al. 
concluded that the excess absorption effect is likely due 
to an intervening galaxy located at z = 0.524, and per- 
formed a more accurate estimate of the spectral indices 



by introducing the absorber at z = 0.524 in the analysis. 
However, the ambiguity on the X-ray spectral behaviour 
of the source is not completely solved, and the debate on 
the origin of X-rays in AO 0235+16 remains still open. 

At 7-ray energies, the two clear detections by EGRET 
in 1994 and 1998 are not enough to constrain the source 
variability range. Only the 1994 signal is sufficiently strong 
to produce a useful spectrum, with a photon index 7 = 
1.88±0.11 (Mukherjee et aHl997: Hartman, private com- 
munication), even if a similar index 7 = 1.85 ± 0.12 was 
derived from the sum of the observations of all viewing 
periods from Cycle 1 up to Cycle 4 (P1234; see Hartman 
et aL iro^ . 

As already noticed for the light curves of Fig. 1, the 
multiwavelength flux changes appear to be superimposed 
to a possible slowly variable "base level" emission, the 
deepest minimum of which is represented by the shaded 
(yellow) area of Figs. 3 and 4: the low-energy part of it 
was determined by extracting the most reliable histori- 
cal minima from well-sampled radio, infrared, and optical 
light curves, and by fitting the corresponding SED with a 
cubic spline interpolation, while the high-energy part was 
analytically modelled according to the SSC theory, under 
the assumption that electron energies vary from 7 = 1 to 
7 = 10^. 

The long-term behaviour of the observed SED vari- 
ability upon the "base level" was modelled as the result 
of the rotation of the steadily-emitting helical-shaped jet 
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previously described. In particular, the low-energy SED 
evolution is consistent with the radio-optical light curve 
modelling, while for the high-energy part two indepen- 
dent cases were considered to take the observation con- 
straints into account: the soft X-ray emission is, in the 
first case, IC radiation, while in the second case it also 
has a synchrotron contribution, which is found to domi- 
nate for small 20°) rotation angles. Difierent-style (and 
colour) curves in Figs. 3 and 4 show the model predictions 
for the SED time evolution in the two different hypothesis, 
as described in more details below. 

Table 1. Helical-jet model input parameters relevant to 
light curve and SED evolution modelling of AO 0235-fl6. 
Col. A: hght curves of Fig. 1 (right panels) and Fig. 5a, and 
SED time evolution of Fig. 3. Col. B: SED time evolution 
of Fig. 4. Values differing in the two models are written 
in boldface. 



Jet geometry 


A 


B 


c 


20° 


20° 




15° 


15° 


a 


360° 


360° 


Jet physics 


A 


B 


r 


10 


10 


7min 


1 


1 


70 


1.995 X 10^ 


6.310 X 10^ 


Cj 








Synchrotron radiation 


A 


B 


ao 


0.5 


0.5 




10^* 


10" 


^^5,0,11 


1.259 X 10^* 


2.512 X 10^" 


^brcak 


0.017 


0.017 


^1,1 ~ h,u 




10-^0 




10"^ 


10"^ 


Cl,I = C2,I 


0.63 


0.63 


Cl,II = C2,n 


9.45 


9.45 


L 


0.1 


0.1 


Ca 


7 


7 




5.012 X 10* 


5.012 X 10* 


IC radiation 


A 


B 


lo 


0.1 


0.1 


Cc 


7 


7 




5.012 


5.012 



The time evolution of the low-energy branch of the 
spectrum, mainly resulting from synchrotron emission of 
the leptonic jet population, is common to both cases. As 
shown in Figs. 3 and 4, the kinematic evolution of the heli- 
cal jet described in Sect. l3.ll can well account for the SED 
changes corresponding to the radio-optical light curve pe- 
riodic behaviour (Fig. 1): as the (p angle increases with the 
rotation and the first, strongly emitting, portion of the jet 
becomes more and more misaligned, the synchrotron peak 
slowly moves towards lower energies, finally plunging into 
the "base level" SED. 

The high-energy part of the SED was studied by tak- 
ing the historical X-ray spectral properties of the source 
into account. As previously mentioned, the shape of the 



X-ray spectrum is extremely variable, exhibiting alterna- 
tively hard and steep indices, and thus raising the con- 
troversial issue of the origin of the source X-ray emission. 
According to the value of a, it could either be synchrotron 
radiation produced by high-energy electrons in the jet, or 
come from IC scattering of soft synchrotron photons (SSC) 
or external radiation (EC) off jet electrons, or be the result 
of a composition of synchrotron and IC effects. 

We hence simulated the high-energy SED time evolu- 
tion in two different cases. 

In the first one (Fig. 3), we assumed that the whole 
X-ray jet emission comes from IC scattering of softer syn- 
chrotron photons off jet energetic electrons [70 w 2 x 10^, 
see Eq. (8)]. 

The X-ray flux variability and the constancy of the 
relevant spectral slope, exhibited by the modelled emis- 
sion as the rotation goes on, reflect the radio SED evolu- 
tion, while the 7-ray SED fluctuations correspond to the 
infrared-optical ones. The parameter setting is shown in 
Col. A of Table 1. 

In the second case (Fig. 4), we supposed the existence 
of a synchrotron contribution to the soft X-ray emission 
to justify the steepest spectra, like those observed by the 
Einstein Observatory in August 1980 and by ROSAT in 
July-August 1993; on the contrary, hard X-rays were still 
assumed to come from IC scattering of softer photons. 
The modelled jet is now characterized by a production of 
synchrotron radiation up to frequencies much higher than 
in the previous case (~ 10"'^^ Hz) in the very first portion 
of the emitting jet. The consequence is the appearance of 
a synchrotron "bump" in the SED of Fig. 4 in the soft 
X-ray band, well visible as long as the jet apex region is 
oriented at small angles {(f> ^ 20°) with respect to the line 
of sight. Moreover, the soft X-ray part of the SED becomes 
much more variable in both flux and slope, thus matching 
better the observed spectra. The 7-ray modelled emission 
also exhibits remarkable new features. It is the result of IC 
processing of a much wider range of frequencies and seems 
to better match observations with relatively low electron 
energies (70 « 6 x lO'^). The flrst SED peak occurs in 
the MeV band and corresponds to the infrared one, while 
the second peak lies in the GeV region and reflects the 
synchrotron X-ray bump, which thus would be the source 
of the seed photons for the EGRET detections. A complete 
description of the model parameters is given in Col. B of 
Table 1. 

4. Modelling the non-periodic events: the 8 GHz 
light curve 

With the aim of investigating the variability of the source 
emission in deeper detail, we attempted a more accurate 
comparison between model and observations for the best- 
sampled radio light curve at 8 GHz. 

The model predicts that, as the helix rotates, the view- 
ing angle of the 8 GHz jet emitting region varies periodi- 
cally, as shown in Fig. 56, implying the light-curve periodic 
behaviour displayed in panel Fig. 5 a. 




Fig. 5. a: Helical-jet modelling of the 8 GHz light curve, b: Time evolution, with the helix rotation, of the viewing 
angle 9 of the jet region emitting radiation observed at 8 GHz. Minima of the viewing angle 9 correspond to maxima 
of the flux, as marked by the vertical dashed line. Only jet emission has been considered here. 
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Fig. 6. a: Cubic spline interpolation through the binned (50 days/bin) 8 GHz light curve (solid line) and "base level" 
flux (dotted line), b: Model time evolution of the viewing angle of the 8 GHz jet emitting region, derived from the 
light curve of panel a in the hypothesis that minor, non-periodic outbursts are due to some distortions of the helical 
structure. 



In order to make the comparison between model and 
data easier, we removed the short-term flaring from the 
observed curve of Fig. 1 c by deriving a binned light curve 
with a bin size of 50 days, and by fitting it with a cu- 
bic spline interpolation (Press et al. I1992|l . as shown in 
Fig. 6a. The minor, non-periodic outbursts interspersed 
among the periodic major ones are clearly visible in the 
curve. We can formulate two hypotheses about their ori- 
gin: (1) they are the effect of some geometric distortions of 
the helical path; (2) they originate from other phenomena 
contributing to the source emission. 



By looking at the interpolation curve in Fig. 6 a, one 
can see that flux minima are less deep than simulated ones 
(Fig. 5a), and that they seem to follow a linear decreas- 
ing trend down to the lowest state, as previously noticed 
for the radio light curves of panels 6, c, d, and g of Fig. 
1. Hence we defined a "base level" for the source emis- 
sion, represented by a straight line from the first datum 
of 1974 (JD = 2442278.6469) to the 1996 deepest mini- 
mum (JD = 2450454.7562), and by a constant value in 
the remaining part of the curve (dotted line in Fig. 6a). 
We then subtracted this trend from the spline, under the 
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assumption that the rcsuhing "cleaned" curve represents 
the jet contribution alone. 

Following hypothesis (1), i.e. supposing the minor out- 
bursts as due to distortions of the helical path, we derived 
the time evolution of the local viewing angle 9 correspond- 
ing to the cleaned curve. The difference between the re- 
sult, displayed in Fig. 66, and the theoretical trend (Fig. 
56) would represent the displacement of the jet path from 
a perfect helix in the jet region of interest. The presence of 
the minor outbursts, missing in the modelled light curve, 
which keeps the mean flux high even when the base level 
is subtracted, is the reason why the average local viewing 
angle 9 is low. 

Under hypothesis (2), the minor outbursts come from 
other phenomena in the source. They might be explained 
as intrinsic variations of the emission, or as due to some 
more complex geometry of the jet structure implying other 
contributions to the flux. By removing the non-periodic 
outburst component from the curve in Fig. 6a, we ob- 
tained the shaded light curve of Fig. 7a, and hence the 
time evolution of 9 shown in Fig. 76, which looks very 
similar to that provided by the model, suggesting that 
the observed periodic behaviour can indeed be accounted 
for by a rotating jet structure fairly stable in time. 

In this latter scenario, we tested the hypothesis that 
the non-periodic features of the light curve originate from 
some MHD instabilities occurring in the jet (e.g. pinching 
modes, which lead to the formation of jet regions char- 
acterized by higher-density magnetic field and stronger 
emissivity), or from other phenomena implying enhanced- 
emission regions (e.g. shock or magnetic reconnection 
events). 

As the helix rotation goes on, these stronger emitting 
regions approach the line of sight, giving rise to outbursts 
in the observed light curve. The presence of enhanced- 
emission zones along the jet was simulated by modifying 
the emissivity fading factor represented in Fig. 2 (dashed. 



red line): Gaussian profiles were added to the curve to 
reproduce, as an example, the main light curve events ob- 
served at 8 GHz over an 8-year time interval centred on 
the 1992-93 periodic outburst. The sequence of zones was 
supposed to be partially evolving during the selected pe- 
riod: the fading factor follows the trend displayed in Fig. 
8c before the outburst, and that in Fig. 8(i after the out- 
burst. The corresponding modelled light curve, shown in 
Fig. 86, well reproduces the observed source behaviour, 
displayed in Fig. 8 a. 

5. Discussion and conclusions 

We found that the helical-jet model is able to describe 
the long-term behaviour of the multiwavelength emission 
of the BL Lacertae object AO 0235+16. Both the peri- 
odic occurrence and the mean shape of the main radio 
and optical outbursts, as well as the corresponding re- 
markable variations of the broad-band SED, can be ex- 
plained in terms of the orientation change of an inhomo- 
geneous, steadily-emitting, rotating helical jet. The radi- 
ation from the different-frequency emitting regions of the 
jet is affected by relativistic beaming, whose amount de- 
pends on the angle between the velocity vector of the emit- 
ting plasma and the line of sight, which changes along the 
helical path and also varies with time. The assumption 
of flow instabilities (or other phenomena implying local 
emission enhancement) in the jet provides a viable inter- 
pretation for the non-periodic outbursts observed in the 
radio light curves. 

The twisting of the jet in a helical structure and its ro- 
tation can originate from the orbital motion of the parent 
black hole in a BBHS, the main signature of which would 
indeed be the periodicity of the source light curves. 

The periodicity of AO 0235+16 is about 2.9 years in 
the host galaxy rest reference frame, taking into account 
the P = Pobs/(l + z) relation with z = 0.94. According 
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Fig. 8. a: Observed 8 GHz light curve during an 8-year period centred on the 1992-93 outburst, b: Modelled 8 GHz light 
curve, under the assumption that non-periodic outbursts originate from the presence of enhanced-emission regions, c: 
Emissivity fading factor before the 1992-93 periodic outburst, d: Emissivity fading factor after the 1992-93 outburst. 



to Begelman et al. H198()|l , this period (assumed to be the 
orbital period) enables to estimate the mass of the primary 
black hole for any given value of the mass ratio M/m 
between the primary and secondary components: M 
Py/^(M/m)3/5 10^ M0. For M/m - 1-100 one can infer 
M ^ 5 X 10^ - 9 X 10^ Mq] the binary separation would be 
in the range 2 x 10"'' - 5 x lO"'' pc. If the jet is emitted by 
the primary, its orbital radius could vary from 5 x 10~^ 
up to lO"'' pc. 

Several VLBI images of the source were produced over 
a wide range of observing radio wavelengths in the past 
years. A large fraction of these maps shows no evidence 
of extended structures apart from the compact core, re- 
gardless of the different resolution (Gabuzda et al. 119921 
Gabuzda & Cawthorne il996) : on the other hand, some of 
them reveal a faint jet north of the core (Jones et al. 119841 
Chu et al.|1996' Shen et al. ll997|l . The presence of a weak 
extension is confirmed by the sub-milliarcsecond maps ob- 



tained with the VSOP at 5 GHz (Frey et al. IMmHi and, 
more recently, with the VLBI at 43 GHz (Jorstad et al. 
I2()()l(l . the latter also displaying a couple of components 
superluminally moving along bent trajectories. 

Our model interprets the jet radio knots as the jet 
regions where the helical pattern presents the minimum 
viewing angle, with a maximization of the beaming effect. 
According to the above orbital radius estimate and tak- 
ing the pitch angle C and the viewing angle of the helix 
axis ^ into account, the observed separation between radio 
knots should be less than 10 /las, well below the resolution 
of the most detailed available maps. All this under the as- 
sumption that the helix pitch docs not vary when moving 
from the "one-turn" emitting region we considered, which 
is likely not true: the helix pitch could be smaller close 
to the black hole and growing outside, up to the observed 
sub-mas scales. 
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A huge observing effort is currently ongoing on this 
source, in order to closely follow its variability be- 
haviour around the time of the next predicted out- 
burst (first half of 2004; Raiteri et al. I^IMIl : op- 
tical and radio telescopes of the WEBT collabora- 
tion (http://www.to.astro.it/blazars/webt/: Villata 
et al.l 2000[ I2002|l are intensively monitoring it since sum- 
mer 2003, together with the Effelsberg 100 m radio tele- 
scope and VLB A. Moreover, the ground-based observing 
effort will be intensified in 2004, during the optical/UV/X- 
ray pointings of the source by the instruments onboard 
XMM-Newton, which will provide more details on the 
shape of the SED in the UV-X-ray band, possibly sheding 
light on the origin of X-rays in AO 0235-1-16. 
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